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Two-dimensional photonic crystal membranes 
provide a versatile planar architecture for inte¬ 
grated photonics to control the propagation of 
light on a chip employing high quality optical 
cavities, waveguides, beamsplitters or dispersive 
elements [l|. When combined with highly non¬ 
linear quantum emitters, quantum photonic net¬ 
works mm operating at the single photon level 
[4] come within reach. Towards large-scale quan¬ 
tum photonic networks IS ED, selective dynamic 
control of individual components and determin¬ 
istic interactions between different constituents 
are of paramount importance [7]. This indeed 
calls for switching speeds ultimately on the sys¬ 
tem’s native timescales [5]. For example, ma¬ 
nipulation via electric fields or all-optical means 
have been employed for switching in nanopho¬ 
tonic circuits [fil 10] and cavity quantum electro¬ 
dynamics studies [TTIIT2]. Here, we demonstrate 
dynamic control of the coherent interaction be¬ 
tween two coupled photonic crystal nanocavities 
forming a photonic molecule fl3l Ml- By using 
an electrically generated radio frequency surface 
acoustic wave we achieve optomechanical tuning 
Ha, demonstrate operating speeds more than 
three orders of magnitude faster than resonant 
mechanical approaches Ha- Moreover, the tun¬ 
ing range is large enough to compensate for the 
inherent fabrication-related cavity mode detun¬ 
ing. Our findings open a route towards nanome- 
chanically gated protocols |17L 18], which hith¬ 
erto have inhibited the realization in all-optical 

schemes mm- 

The on-chip photonic molecule (PM) studied here con¬ 
sists of two L3-type missing hole cavities, labeled C\ and 
C 2 , defined in a two-dimensional photonic crystal mem¬ 
brane. Its static coupling strength J is known to depend 
exponentially on the separation between the two cavi¬ 
ties El- In our sample the cavities are offset symmet¬ 
rically by d = 5 lattice constants (a = 260 nm) along 
each primitive direction of the photonic crystal (PC) lat¬ 
tice (see Supplementary information), as can be seen in 


the scanning electron micrograph in Fig. [l] (a). The fi¬ 
nite coupling strength J leads to the formation of two 
normal modes; a bonding mode M_ and an antibond¬ 
ing mode M_|_. Note that the mode indices in M± refer 
to the respective normal mode frequencies uj± and, thus, 
are opposite to the spatial symmetry of the mode. In the 
most general case, the two cavities C\ and C 2 exhibit res¬ 
onance frequencies uq and 002 , split by a finite detuning 
A = U 02 — cq. Thus, the normal mode frequencies can be 
expressed by 

UJ ± = w 0 ± i a/a 2 + J 2 (1) 

with ujq = (uq + 0^2) / 2 being the centre frequency. In 
Fig. [l] (b) we present two emission spectra of embed¬ 
ded quantum dots (see methods) recorded for spatially 
exciting either C\ (blue) or C 2 (red) of a typical PM, de¬ 
noted as PM1. For PM1, each cavity exhibits a distinct 
single mode, split by co± = ujq db |y/Ag + J 2 , with A 0 
being the static detuning. Therefore, these as-fabricated 
and nominally identical cavities are efficiently decoupled 
due to inevitable imperfections during nanofabrication. 
To overcome this fundamental limitation and achieve dy¬ 
namic control of the individual cavity resonances we em¬ 
ploy an optomechanical approach using radio frequency 
(rf) surface acoustic waves (SAWs) [15l : '2J]. As indicated 
in the schematic illustration of Fig. [l] (a), these quasi- 
monochromatic acoustic phonons are generated by an in¬ 
terdigital transducer (IDT) (see methods), and propagate 
along the x-axis of the PM and dynamically deform the 
individual cavities. The amplitude of this mechanism, 
A/2tt > 250GHz, exceeds both the static detuning, Aq, 
and the coupling strength, J, which enables us to dynam¬ 
ically tune the PM completely in and out of resonance. 
For our chosen geometric arrangement, the cavities are 
separated by a distance Ax along the SAW propagation 
direction, the modulations of uq and uq are phase shifted 
by 4> 12 = Ax ■ wsaw/csaw, with w S aw and c S aw be- 
ing the SAW angular frequency and speed of sound, re¬ 
spectively. We deliberately set the acoustic wavelength, 
Asaw, such that it is commensurate with the spatial sep¬ 
aration of the two cavity centers 2Ax = Asaw- As il¬ 
lustrated in the schematic in Fig. [I] (c), for the selected 
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FIG. 1. Photonic molecule with acousto-mechanical tuning mechanism, (a) SEM image of the photonic molecule and schematic 
of the device layout. The IDT launches a SAW which propagates across the photonic molecule. Scale bar, 2 pm (IDT not 
to scale), (b) PL spectra showing the two fundamental modes of the photonic molecule and optical microscope images of the 
excitation laser spot. The blue and the red curve were obtained by exciting on cavity C\ and C 2 , respectively. Due to a finite 
static detuning Ao > J, the modes are confined within their respective cavities. The center frequency ujq is 2n x 305 THz. The 
cavities can be spatially resolved in our optical setup, allowing selective excitation of the individual cavities. Scale bar, 5 pm. 
(c) Acousto-mechanical tuning. The resonance frequencies of the cavities are modulated by the SAW induced strain field 
(amplitude exaggerated). Setting half the SAW wavelength equal to the cavity separation, the cavities can be tuned relative 
to each other, giving rise to a time dependent detuning A (t) — Ao + A mo d sin^sAwt)- (d) Evolution of the modes over 
one SAW period. The dashed lines represent the single cavity modes Mi and M 2 , and the solid lines the normal modes M+ 
and M-. The color scale corresponds to the decomposition into single cavity modes, with green denoting the symmetric and 
antisymmetric superposition on resonance. The static detuning Ao and the coupling strength J are intrinsic parameters of the 
photonic molecule. The tuning amplitudes A i j2 and the phase shift 0i2 are determined by the amplitude and wavelength of the 
SAW, respectively, (e) FDTD simulations of the mode profiles ( E y ). For the case of detuned cavities (A = 3J), the modes are 
localized within the individual cavities. On resonance (A = 0), one obtains the symmetric (bonding) and the antisymmetric 
(antibonding) mode. 


phase, the maximum (minimum) of the SAW expands C\ 
(compresses C 2 ), which in turn gives rise to a red (blue) 
shift of their respective resonances [15 . The propagation 
of the sound wave along the x-axis of the PM leads to a 
time-dependent detuning of the two cavities 

A(£) = A 0 + A mo d sin(cdsAW^)- (2) 

The amplitude of this modulation, A mo d = 
y/A\ + A%- 2A^A~ 2 cos( 012 ), depends on the amplitudes 
of the modulations of the individual cavities, A\ and A 2 , 


and is maximum for 12 = 180°. For A mo d > A 0 , the 
detuning passes through zero, which results according 
to equation 0 in an avoided crossing of the normal 
modes. In Fig. |T] (d), the resulting temporal evolutions 
of the normalized mode frequencies (uj — cjo/J) obtained 
using equations 0 and 0 are evaluated over one 
acoustic cycle of Tsaw for a typical set of experimentally 
achievable parameters (J = 1.2 • Ao, A\ = 1.5 • Ao, 
A 2 = 1.9 • Aq and </>i 2 = 162°). The dashed lines show 
the time-evolution of the non-interacting (J = 0) modes 
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Mi (blue) and M 2 (red). In strong contrast, the normal 
modes M + and M_ start with initially Mi-like (blue) 
and M 2 -like (red) single-cavity character and develop 
to fully mixed symmetric (bonding) and antisymmetric 
(antibonding) character (green) at resonance. At this 
point, the coupling strength J of the PM can be deduced 
directly from the splitting of the avoided crossing. After 
traversing through the avoided crossing, the initial 
character of the modes is exchanged, with M + and 
M_ possessing M 2 -like (blue) and Mi-like character, 
respectively. Over the duration of one full acoustic cycle 
the two modes are brought into resonance at two distinct 
times, giving rise to two avoided crossings, restoring the 
initial configuration. We performed full finite difference 
time-domain (FDTD) simulations to confirm and in 
particular quantify the dynamic coupling. We evaluate 
the calculated profiles of the electric field component 
in the plane of the photonic crystal membrane and 
perpendicular to the SAW propagation (E y ) in Fig. [l] 
(e). In the two left panels, the detuning is set to three 
times the coupling strength (A = 3J) and, thus, the two 
modes remain well localized in one of the two cavities. 
For the resonance case, A = 0, shown in the two right 
panels of Fig. [I] (e), the E y exhibit the characteristic 
symmetric and antisymmetric superpositions for the M_ 
and M + modes, respectively. For our sample geometry 
the FDTD simulations predict J s i m /27r = 170 GHz. 
This value is set by the chosen separation, i.e. barrier 
thickness, between the two cavities. 

To experimentally verify the dynamic control of a PM 
by a cjsaw/27t = 800 MHz SAW we performed strobo¬ 
scopic spectroscopy which ensures optical excitation at 
a well-defined time during the acoustic cycle. The emis¬ 
sion of the PM was analysed in the time-domain and as a 
function of relative emission frequency uj — ujq . We mea¬ 
sured time dependent emission spectra of the two cavities 
of PM1 for three characteristic modulation amplitudes. 
For a weak modulation, A mod < A 0 , shown in Fig. [ 2 ] 
(a), C\ and C 2 show the expected phase-shifted sinu¬ 
soidal spectral modulations centred (dashed horizontal 
lines) at their unperturbed resonances which decay with 
a characteristic time constant of ~ 1.3 ns. As the mod¬ 
ulation is increased to A mo d = Aq (Fig. [ 2 ] (b)) the two 
single cavity modes are brought into resonance at one dis¬ 
tinct time to = 0.6 ns during the acoustic cycle. At this 
time, coherent coupling leads to the formation of bonding 
and antibonding normal modes. This directly manifests 
itself in the experimental data due to the emergence of 
new emission features stemming from the spatial delo¬ 
calization of the bonding and antibonding modes. For 
the initially lower frequency cavity (Ci), a new signal 
appears at the frequency of the normal mode M + , which 
was initially confined within the other cavity (C 2 ). The 
initially higher frequency cavity (C 2 ) exhibits precisely 
the opposite behaviour, with the normal mode M_ ap¬ 
pearing at time £q. The normal modes are split by the 


coupling strength J/2tt = 164 ± 10 GHz, very close to 
the value expected from our FDTD simulations. For fur¬ 
ther increased detuning A mo d > Ao, the two modes are 
brought into resonance at two distinct times, ti and £ 2 , 
during the acoustic cycle. After the first resonance at 
£ 1 , coupling is suppressed, A > J, both modes are ef¬ 
fectively decoupled and their single cavity characters are 
exchanged compared to the initial configuration. The 
lower frequency mode M— , which is initially Ci-like, is 
switched to CVlike character after £ 1 , and vice versa. At 
the second resonance at £ 2 , the system is reverted to its 
original configuration at the beginning of the acoustic 
cycle. This sequence of two time-offset coupling events 
gives rise to the experimentally observed anticorrelated 
time-evolutions of the emission of C\ and C 2 . We ex¬ 
tracted the time-dependent emission frequencies detected 
from the two cavities of PM 1 and plotted them (C\ cir¬ 
cles, C 2 triangles) for the three modulation amplitudes in 
Fig. [ 2 ] (d-f). All experimental data are well reproduced 
by the simple coupled mode model given by equations 
0 and 0. From the set of measurements with vary¬ 
ing tuning amplitude we obtain mean values for the free 
parameters of the model, (J)pmi/27t = 170 ± 10 GHz, 
(Ai/A. 2 )pmi = 0.75 i 0.05 and (0 i2)pmi = 135 =b 10°. 
Moreover, we note that the observed phase shift unam¬ 
biguously confirms efficient coupling of the SAW into the 
sub- As aw membrane. Its small deviation from the ideal 
value of 180° arises from the finite difference in the phase 
velocity of the acoustic wave within the membrane and 
the region of the transducer. The results from this model 
are plotted as lines and the character of the mode is 
color coded. Clearly, our experimental data are in ex¬ 
cellent agreement with the normal mode model for all 
three modulation amplitudes. Both the time-dependent 
spectral modulation as well as the switching of the char¬ 
acter of the modes are nicely reproduced. 

Such behaviour was experimentally confirmed for 
five different, nominally identical PMs for which we 
evaluated the mean of their respective key parameters. 
The experimentally observed static detuning Ao and 
coupling strength J are summarized in Fig. 1(a)- 
While the coupling strength, J/2tt = 192 ± 26 GHz, 
(blue) does not vary from PM to PM, the values of the 
static detuning, Ao, shows a pronounced scatter ranging 
between —300 GHz and +900 GHz. These observations 
are in fact expected. J exponentially depends on the 
symmetric inter-cavity offset, d, and is thus robust and 
insensitive to small deviations from the ideal geometry 
due to fabrication imperfections. In strong contrast, 
the absolute resonance frequencies of the two cavities 
forming the PM are highly sensitive to these inherent 
and inevitable deviations from the nominal geometry. 
The resulting fluctuations of the cavity resonances 
reflect themselves in the observed pronounced variation 
of Ao- Indeed, we observe pronounced coupling effects 
(|J| > |Aq|) only for one single as-fabricated PM, 
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FIG. 2. Experimental characterization for three different modulation amplitudes, (a-c) Time resolved PL maps, measured 
on cavity C\ (upper panels) and C 2 (lower panels), (d-f) Mode frequencies extracted from the measurements on C\ (circles) 
and C 2 (triangles), fitted with the coupled mode model. For A mo d < Ao (a and d), the modes are effectively decoupled. Their 
frequencies are modulated sinusoidally by the SAW with a phase difference of 0i2 = 160°. For A mo d = Ao (b and e), there 
is a time to at which the cavities are in resonance. At this point the modes are split by the coupling strength J and they are 
delocalized over the photonic molecule, thus, both modes can be observed simultaneously on each cavity. When A mo d > A 0 (c 
and f), the modes exhibit an avoided crossing. This happens twice in each SAW cycle, at times t\ and £ 2 - 


labeled PM5. The corresponding experimental data 
is presented in the Supplementary information. The 
exponential dependence of J as a function of d is 
nicely confirmed by a best fit (line) to values calculated 
by FDTD (symbols) presented in Fig. [ 3 ] (b). The 
experimentally observed distribution of J exp and the 
inter-cavity separation d exp derived from J exp and the 
FDTD simulation are indicated by the shaded horizontal 
and vertical bars. Clearly, the measured J exp and its 
derived d exp matching perfectly the calculated value of 
^sim/27r = 170 GHz and the nominally set d = 5a. These 
narrow distributions centred around the calculated and 
nominal values are expected since d is large compared 
to typical imperfections in the nanofabrication. 


pled nanophotonic system. When combined with opti¬ 
cal non-linearities our tunable PM paves the way to dy¬ 
namically controlled high-fidelity entanglement genera¬ 
tion m and distribution on a chip |25( [26] . For larger 
switching rates as required for Landau-Zener-transit ion 
based gates PH, the underlying optomechanical coupling 
could be enhanced further by direct SAW excitation of 
localized vibronic modes of PCM nanocavities [23] or 
shaped SAW waveforms [T8] , Finally we note, that our 
approach could be directly scaled up to large arrays of 
coupled cavities EH El] or employed in superconduct¬ 
ing two-level systems, which have recently been strongly 
coupled to single SAW quanta [29 . 


In summary, we demonstrated dynamic optomechani¬ 
cal control of coherent interactions in a prototypical cou- 
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FIG. 3. Statistical variation of PM properties and comparison 
with FDTD simulation, (a) Coupling strength J and static 
detuning Ao measured on five nominally identical PMs. (b) 
FDTD simulation of the coupling strength for different cav¬ 
ity separations d. The experimentally determined coupling 
strength is accurately reproduced by the simulation. 


METHODS 
Sample structure. 

We start by fabricating the photonic crystal mem¬ 
branes from a semiconductor heterostructure grown by 
molecular beam epitaxy. This heterostructure consists of 
a 170 nm GaAs layer with self-assembled InGaAs quan¬ 
tum dots (QDs) at its centre, on top of a 725 nm thick 
Alo. 8 Gao. 2 As sacrificial layer. The PM structure is de¬ 
fined by electron beam lithography and transferred into 
the heterostructure by ICP-RIE etching. In a wet chem¬ 
ical etching step using hydrofluoric acid we removed the 
sacrificial layer to release a fully suspended membrane. 
The PMs are deliberately designed to be off-resonant 
with the QD emission to achieve a sufficiently long de¬ 
cay time of the cavity emission. Thus, we can observe 
spectral modulations over more than one acoustic cy¬ 
cle in the time domain, which would be impeded by 
Purcell-enhanced resonant QD exciton recombinations 
m- IDTs were defined using electron beam lithogra¬ 
phy and metallized with 5 nm Ti and 50 nm A1 in a 
lift-off process. The finger period is 3.83 pm, resulting in 
a resonance frequency of 800 MHz at 5 K. 


Optical spectroscopy. 

For stroboscopic photoluminescence spectroscopy of 
the PM, off-resonant QDs are excited by an externally 
triggered diode laser emitting 90 ps pulses at a wave¬ 
length of 850 nm which is focused to a 1.5 pm spot by 
a NIR 50 x microscope objective. The emission is dis¬ 
persed by a 0.75 m imaging grating monochromator. A 
liquid N 2 -cooled Silicon charge coupled device and a fast 
(< 50 ps rise time) Si-single photon avalanche detector 
(SPAD) are used for time-integrated multi-channel or 
time-resolved single channel detection, respectively. The 


sample is cooled to T = 5K in a Helium-flow cryostat 
with custom built integrated rf connections. 


Coupled mode model. 


We treat the PM as a model system of two coupled 
cavities. The single cavity modes M\^ have the complex 
amplitudes < 2^2 and the frequencies uj i ? 2 - In the presence 
of finite coupling and the absence of dissipation, the time 
evolution is given by 


da\ . .J 

—— = lujxai + i-a 2 

dt 2 

da 2 . .J 

—— = iu 2 a 2 + i — a\ 
dt 2 


with a real coupling constant J. The resulting nor¬ 
mal modes have the frequencies uj± = cjo d= | a/A 2 + J 2 
[Equation ([!])] with centre frequency uj 0 = ^(uji+uj 2 ) and 
detuning A = uj 2 — cji . We note that since J > 0 the low 
frequency mode uj- and the high frequency mode uj+ cor¬ 
respond to symmetric and antisymmetric superpositions 
of the uncoupled, single cavity modes. 
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